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We investigate the scaling properties of products of the exponential of birth-death processes 
with certain given marginal discrete distributions and covariance structures. The conditions on 
the mean, variance and covariance functions of the resulting cumulative processes are interpreted 
in terms of the moment generating functions. We provide four illustrative examples of Poisson, 
Pascal, binomial and hypergeometric distributions. We establish the corresponding log-Poisson, 
log-Pascal, log-binomial and log-hypergeometric scenarios for the limiting processes, including 
their Renyi functions and dependence properties. 
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1. Introduction 

Fully developed turbulence has been characterized by certain universal properties such 
as the scaling behavior 

across a distance I of the gth-order moment of velocity fluctuations, or 

E(e,)«^r(«' 

for the qth-order moment of locally averaged energy dissipation e; over a ball of size I. 
Kolmogorov's refined similarity hypothesis [20] leads to the relationship (Frisch [11]) 

Ciq) = l+riq/3). 
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2009, Vol. 15, No. 2, 508-531. This reprint differs from the original in pagination and 
typographic detail. 



1350-7265 © 2009 ISI/BS 



Multifractal scaling of products of birth-death processes 



509 



She and Leveque [30] proposed a scaling model which predicts 



( 



2 



) 



C(9) = ^ + 2 1 



3 



DubruUc [8] explored further properties of the model and showed that its probability 
density function is related to the log-Poisson statistics of local non-dimensional energy 
dissipation. She and Waymire [31] proposed that the statistics of the inertial range of 
fully developed turbulence can be described by random multiplicative cascades, showed 
that these random cascades arc log-Poisson and re-established the scaling model of She 
and Leveque [30]. 

This paper will provide a construction of the log-Poisson random cascade and certain 
other models of the same type, such as log-Pascal, log-binomial and log-hypergeometric 
cascades, in the framework of multifractal products of birth-death processes. Multifractal 
products of stochastic processes were defined in Kahane [16, 17] and further investigated 
in Mannersalo et al. [22]. In the present paper, products of independent rescaled copies 
of a "mother" process are considered where the mother process is in the form of the 
exponential of a birth-death process with certain given marginal discrete distribution 
and covariancc structure. The birth-death process is an important example of a homo- 
geneous Markov chain {X(t), t > 0} on the state space S — {—1,0, 1, . . .}. Here X{t) 
can be considered as the size of a population at time t, which fluctuates according to 
the following rule: If at time t the chain is in a state i € S, in one transition it can go 
only to i — 1 or I + 1 . The transition from state i to state z -I- 1 indicates a "birth" in 
the population, whereas the transition from i to i — 1 indicates a "death" . A classifica- 
tion of birth-death processes was given in Karlin and McGregor [18]. A general theory 
of birth-death processes can be found in Dynkin and Yuskevich [9], Karlin and Taylor 
[19] and Bhattacharya and Waymire [6]. For the cases of Poisson, Pascal, binomial and 
hypergeometric distributions of these processes, we construct log-Poisson, log-Pascal, 
log-binomial and log-hypergeometric scenarios with non-trivial singularity spectra. Our 
method is based on an application of proper estimates of the maximal increments of the 
process in conjunction with the orthogonal expansion of its transition distribution. 

There are many constructions of random multiplicative cascades ranging from the 
simple binomial cascade to measures generated by branching processes and the compound 
Poisson process (see Kahane [16, 17], Molchan [23], Falconer [10], Barral and Mandelbrot 
[5], Riedi [28], Morters and Shieh [24, 25, 26] and Shieh and Taylor [32]). Many of these 
multifractal models were not designed to cover dependence structure and a natural form 
of the singularity spectrum (see Novikov [27] and Riedi [28], e.g.). Dependent cascades 
were first considered in Waymire and Williams [34], [35], including a general theory for 
Markov dependent cascade generators. Jaffard [15] showed that Levy processes (except 
Brownian motion and Poisson process) are multifractal; but since the increments of a 
Levy process are independent, this class excludes the effects of dependence. Moreover, 
Levy processes have a linear singularity spectrum while real data often exhibit a strictly 
concave spectrum. In this paper, we pay attention to the strong correlation and nonlinear 
form of the singularity spectrum of a class of random cascades. 
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Our exposition relies on some results of Mannersalo et al. [22] on the basic properties 
of multifractal products of stochastic processes, but we provide a new interpretation of 
the conditions on the mean, variance and covariance functions of the resulting cumulative 
processes in terms of the moment generating functions. This interpretation is more useful 
for our development. We describe the behavior of the gth-order moments and Renyi func- 
tions, which are nonlinear, hence displaying the multifractality of the limiting cumulative 
processes. A property on their dependence structure, leading to their possible long-range 
dependence, is also obtained. The development of Mannersalo et al. [22] can be applied 
to a large class of stationary processes such as stationary diffusion processes. As applica- 
tions, Mannersalo et al. [22] looked at continuous-time two-state Markov processes, while 
in the present paper we also consider infinite state space and some new scenarios such 
as those cited above. 



2. Multifractal products of stochastic processes 

This section recaptures some basic results on multifractal products of stochastic pro- 
cesses as developed in Kahanc [16, 17] and Mannersalo et al. [22]. We provide a new 
interpretation of their conditions based on the moment generating functions, which is 
useful for our exposition. 

We introduce the following conditions: 

A'. LetAW(i), te[0,l], i = l,2,...,be a sequence of independent, strictly station- 
ary, positive stochastic processes such that, for all t, ti, t2 S [0, 1] and i = 0, 1, 2, . . . , the 
following assumptions hold: 

EA«(t) = l, (2.1) 

VarA(')(t) =cri <oo, (2.2) 

Cov(AW(ti),AW(t2))-i?A(ti-t2) = aip,(ti-t2), P.(0) = 1. (2.3) 

We consider the following setting: 

A". Let A^''' be independent rescalcd copies of some measurable separable mother 
process A, that is, 

A['^(t)=A«(t6'), te[0,l], 4 = 0,1,2,..., 

where the scaling parameter b > 1, EA(i) = 1, and = denotes equality in finite- 
dimensional distributions. 

Moreover, in the examples of Section 4, the stationary mother process satisfies the 
following conditions: 

A'". For t G [0,1], let A{t) = exp{X{t)}, where X{t) is a stationary process with 
EX^{t)<oo, 

Cov(X(ti),X(t2)) = Rxih - t2) = aj,rx{ti - ^2), rx{0) = 1. 
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We assume that there exist a marginal probabiUty mass function pe{x) and a bivariate 
probabihty mass function pg(xi,X2',ti — t2) such that the moment generating function 

M(C)=Eexp{CX(i)} 
and the bivariate moment generating function 

M(Cl,C2;il-t2)=Ecxp{ClX(ti)+C2X(t2)} 

exist. Here, 9 is the parameter vector of the mass function of the process X{t). 
Under the conditions A'-A'", the assumptions (2.1)-(2.3) can be rewritten as 

EA«(t) = M(l) = l; 
Var A['^(t) M{2) - 1 cri < oo; 
Cov(Ai')(ti), a['' {t2)) = 1; (ti - t2)h') - 1, 6 > 1. 
We define the finite product processes 



A„(t) = nA«(t)=expK]x«(i6^ 



i=0 k i=0 ) 

and the cumulative processes 

A„(t)= / A„(s)ds, n = 0,l,2,..., 

where X«(i),i = 0,... , n, . . . , are independent copies of a stationary process i > 0. 

We also consider the corresponding positive random measures defined on Borel sets B 
of [0,1]: 



^JLn{B) = / A„(s)ds, n = 0,l,2,, 

JB 



Kahane [17] proved that the sequence of random measures ^„ converges weakly almost 
surely to a random measure Moreover, given a finite or countable family of Borel 
sets Bj on [0,1], it holds that lim„_,oo /x„(i3j) = fi{Bj) for all j with probability one. 
The almost sure convergence of An{t) in countably many points of [0, 1] can be extended 
to all points in [0, 1] if the limit process A(t) is almost surely continuous. In this case, 
lim„^oo^„(i) = A{t) with probability one for all i € [0,1]. As noted in Kahane [17], 
there are two extreme cases: (i) An{t) A{t) in Li for each given t, in which case 
A{t) is not almost surely zero and is said to be fully active (non-degcncratc) on [0,1]; 
(ii) An{l) converges to almost surely, in which case A{t) is said to be degenerate on 
[0, 1]. Sufficient conditions for non-degeneracy and degeneracy in a general situation and 
relevant examples are provided in Kahane [17] (equations (18) and (19), respectively.) 
The condition for complete degeneracy is detailed in Theorem 3 of [17]. 
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The Renyi function, also known as the deterministic partition function, is defined as 



ioel4"'l 



= liminf(-i)log2E^M''(4"'), 



k=0 

where /^"^ = [fc2-", (/c + 1)2""], A: = 0, 1, . . . , 2" - 1, is its length and logf, is log 

to the base b. 

Remark 1. The multifractal formalism for random cascades can be stated in terms of 
the Legendrc transform of the Renyi function: 

T*(a) =min((7a — T(q)). 
In fact, let /(a) be the Hausdorff dimension of the set 

= <j < e [0, 1] : hm^ — — ^ "^^^^ = 



k 



where l["\t) is a sequence of intervals ij^^ that contain t. The function /(a) is known 
as the singularity spectrum of the measure /i, and we refer to /i as a multifractal measure 
if /(a) 7^ for a continuum of a (Lau [21]). In order to determine the function /(a), 
Hentschel and Procaccia [14], Frisch and Parisi [12] and Halsey et al. [13], for example, 
proposed to use the relationship 

f{a)^T*{a). (2.4) 

This relationship may not hold for a given measure (see, e.g., Taylor [33]). When the 
equality (2.4) is established for a measure /i, we say that the multifractal formalism 
holds for this measure. 

Mannersalo et al. [22] presented the conditions for L2-convergence and scaling of mo- 
ments. 

Theorem 1 (Mannersalo, Norms and Riedi [22]). Suppose that the conditions A! -A!" 
hold. 

If, for some positive numbers S and 7, 

exp{-<51rl} < pir) = ^^^^^^^^^^^ < ICr^ , (2-5) 

then An {t) converges in L2 if and only if 

b>\ + al = M(2). 
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If An{t) converges in Li, then the limit process A{t) satisfies the recursion 

A{t)^^j\{s)dA{bs), (2.6) 

where the processes A(t) and A{t) are independent, and the processes A{t) and A{t) have 
identical finite- dimensional distributions. 

If A(t) is non-degenerate, the recursion (2.6) holds, A{1) e Lq for some q > and 
c(g, &~") < oo, where c((7,t) = Esup^gjQ jj |A'^(0) — A''(s)|, then there exist con- 
stants C and C_ such that, for all t G [0, 1], 

^^q-log,EA«(t) <E^<J(^) <^^9-log,EA''(t)^ (2.7) 

which will be written as 

(t) ~ EA' [t) ^ 

If, on the other hand, A{1) G Lq, q> 1, then the Renyi function is given by 

T{q) = q-l- logfcEA«(f) = q-l- log, M{q)- (2.8) 
If A(t) is non- degenerate, A{1) G L2, and A{t) is positively correlated, then 

VarA(<) > Var /" A(s)ds. (2.9) 
Jo 

Hence, if /q A(s)ds is strongly dependent, then A{t) is also strongly dependent. 

Remark 2. The result (2.7) means that the process A{t), tG [0,1] with stationary 
increments behaves as 

logE[^(i + S)- A{t)Y - Kiq) log (5 + Cq (2.10) 

for a wide range of resolutions S with a nonlinear function 

K{q) = q - log,EA''(i) = q - log, M{q), 

where Cq is a constant. In this sense, the stochastic process A{t) is said to be multifractal. 
The function K[q), which contains the scaling parameter b and all the parameters of 
the marginal distribution of the mother process X{t), can be estimated by running the 
regression (2.10) for a range of values of q. For the four examples in Section 4, the explicit 
form of K{(i) is obtained. Hence these parameters can be estimated by minimizing the 
mean square error between the K[g) curve estimated from data and its analytical form 
for a range of values of q. This method has been used for multifractal characterization 
of complete genomes in Anh et al. [2]. 
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3. Geometric birth— death processes 

We consider the onc-dimensional birth-death process X ~ {X{t), t >0} on the state 
space 5 = 1— 1,0,1,...}. This process is defined by its birth-death rates Xi,fj.i,i€ S 
(Karlin and McGregor [18]), that is, X is a Markov process with stationary transition 
probabihties 

p,, (t) = P{xit + s)=j\x{s) = i}, s, 

which are independent of s. In addition, we assume that Pij{t) satisfy 

P^^^+l{h) = X,h + o{h), h^O,ieS; (3.1) 

P,,,_i(/i) = /iJi + o(/i), h^O,ieS; (3.2) 

P,4h) = l-{X,+^i,)h + o{h), h^O,ieS; (3.3) 

P,,,(0) = J^; P_i,_i(t) = l, P_i,,,(i)=0, t>0, ^^-l. 

The order o{h) may depend on i, and thus we will use the notation o{h;i) when it 
becomes necessary. We assume > 0, Ao > 0, Ai, yU^ > 0, i > 1. 

If yUo > 0, then we have an absorbing state —1; once the process enters —1, it can never 
leave it. If /xq = 0, we have a reflecting state 0. After entering 0, the process will always 
go back to state 1 after some time. In this case, the state —1 can never be reached and 
so can be ignored, and we take S = {0, 1,2,.. .}. 

We may then suppose that — 0, so that the state —1 is ignored. We assume that 

oo oo ^ 

71":= V^TTfc <oo, V"- =oo, 

AfcTTfc 

fc=0 k=0 



where the potential coefficients tt^ are given by 

AqAi • • • Ai_i 



« > l,7ro = 1, 



so that the process is ergodic (see Karlin and Magregor [18], Theorem 2). Then, 



limP,,,(i):=p, = ^, j>0, 



exist and are independent of the initial state i. We write the stationary distribution as 
P = (jjj, j = 0, 1,2,...). The infinitesimal generator A of the process is given by 

Afit) = A J(z + 1) - (A, + ^l,)f{^) + M J(i - 1) 

for all bounded real- valued functions /. 

We denote by hi^, P) the Hilbert space of functions f{j),j = 0, 1, . . . , such that 

oo 

U)P J <0O. 
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The spectral analysis of birth-death processes is based on the birth-death polynomials 
Qn{x),n = 0, 1, 2, ... , which are defined by the recursive relations 

- xQn{x) = UnQn-lix) - (A„ + ^i„)Qnix) + XnQn + lix), ?1 = 0, 1, 2, . . . , (3.4) 

with (5_i(x) = 0, Qo{x) — 1. Karlin and McGregor [18] proved that there exists a positive 
Borel measure (j){dx) with total mass 1 and with support on the non-negative real num- 
bers, called the spectral measure of the process, such that the transition probabilities are 
represented in the spectral form as 



The polynomials Qnix) are orthogonal with respect to (j). For certain choices of birth- 
death parameters, wc can apply the theory of classical orthogonal polynomials with 
respect to discrete distributions. In the case /iq = 0, which wc assume in this paper, we 
have Qn{0) = 1, n. > and that is i/i-atomic; the stationary distribution P is then given 



The following proposition plays a key role in our exposition. Wc note that the assump- 
tions (ii) and (iii) in Section 4 of Manncrsalo et al. [22] generally do not hold in the case 
that the state space S is infinite; indeed their assumption (ii) assumes that the jump rate 
v{x) is uniformly bounded from below and from above for all x; this assumption is not 
suitable for an infinite state space, as our present Section 4 shows. We also remark that 
the condition (3.5) below, which is needed for the validity of Proposition 1 for an infinite 
state space, may not be the best, and the possibility of improving it is left for a future 
study (wc thank the referee for noticing this). In the following, denotes expectation 
with respect to the stationary distribution P. 

Proposition 1. We assume, besides the strict stationarity, that the small order of mag- 
nitude in the transition probabilities (3.1)-(3.3) is such that o(h;k) = k^o{h;l) for all 
fc = 1, 2, 3, . . . , for some S > 0. Moreover, we assume that 




= TT. 




byp,=7r,0({O}). 



OO 




(3.5) 



Let g{x) c'^^^ where c G M. Then, for any b>\, 




(3.6) 



Proof. Wc first note that 




(3.7) 



k 
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Since the process is cadlag and integer- valued, we have, for each < > 0, 
E f max \g{X{s)) - 9{X{0})\ \ X{Q) = k 

\0<s<t 

<e( J2 \9{X{s))^giXis-)) X{0)^k 

\<s<t, AX^^O 

where AXg = X{s) — X{s—),0 < s <t, and those s for which AXg ^ are countably 
many. The defining property of birth-death processes imphes that 

P{AXs + 0|X(s-)} < const X (Ax(.-) + Ma-(.s-)), 

where the constant is taken to be the expected jump time +1, which is independent of 
a. We note that, by the definition of .g(-), 

< \g{X{s)) - g(X{s~))\ < e«Siiax{|e« - 1|, |1 - e"«|} 

when X{s-) = h and AX^ ^ 0. As U 0, the joint probabihties of {X(s-) = h,X{0) = k}, 
h varying over S — {fc}, are all of order o(t; fc), except when h = k + I and h = k — 1, in 
which case they are of order \kt and /ifct, respectively. Therefore, for large n and by the 
strong Markov property of the process, which starts afresh at each jump time, we have 

Ef max \g{X{s))-g{X{m\X{0)^k 



<E 



<E 



( \9{X{s))^g{X{s-))\\x{Q)^k\ 

( - g{X{s~))\ I X{s~)) I X{Q) ^ k\ 



<Ex(o)=fc E E(|.9(X(s))-5(X(s-))||x(s-) = fc±l) 

\0<s<h-",AXa5^0 / 

+ o(6-";A:) 

<C,(e«'=(Afe+//fe) + ^')^'""- 

In the above inequalities, we have arranged, without indexing explicitly, the jump times 
in successively increasing order, and have made use of the assumption o(/i; k) — k^o{h; 1) 
for all k = 1,2,3, ... . We note that the expected number of the jumps is at most 1 on 
[0,6~"], and thus the summation of jumps in the above is indeed reduced to one term 
only. In view of the conditional expectation formula (3.7) and the definition of Sq, which 
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we assume to be finite, we get 



max \giXis)) g(X(0))| ] < (C, ■ S,) ^ 6"". 

n=0 ' ' ' n=0 

This completes the proof of the proposition. □ 



Remark 3. We usually consider g > 0, in which case the finitcness of Sq depends on the 
growth of the sequences Afc , ■ We have included here the case q < since it incurs no 
additional difficulty. 



4. Multifractal scenarios 

This section introduces four illustrative multifractal scenarios. The mother process will 
take the form K{t) = exp{X(<) — cx}, where X{t) is a stationary birth-death process and 
Cx is a constant depending on the parameters of its marginal distribution. This form is 
needed for the condition EA(<) = 1 to hold. 



4.1. Log-Poisson scenario 

The log-Poisson statistics in fully developed turbulence were discussed in DubruUe [8] and 
She and Waymire [31]. In this section, we provide a related model, namely the log-Poisson 
scenario for multifractal products of stochastic processes. 
B'. Consider a mother process of the form 

A(t) = exp|x(t)-^(e-l)|, (4.1) 

where X{t)^ t > is a stationary birth-death process with marginal Poisson distribution 
Poiff), A > 0,yu > 0, and rates 

A„ = A, pin= fJ-n, n > 0. 

The covariance function of the process X{t) then takes the form 

Rx{t)^-rx{t), rx(t)=e-'^*. 

Under condition B', we obtain the following moment generating function: 

Af(C) = Eexp|c(^X(i)-^(e-i; 



:exp|-C-(e-l) + -(e?-l)l, C^: 
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It turns out that, in this case 

■og.EA.W^ -'^^'''°-;4V^'"°'-" . ,>0. (4.2, 

We can formulate the following: 
Theorem 2. Suppose that condition B' holds. Then, for any 



the stochastic processes 



6>exp|^(e-l)2|, 
Mt)= / TTA«(s6^")ds, te[0,l], 



3=0 

converge in L2 to the stochastic process A{t), i G [0, 1] as 00 such that, if A{1) G Lq 
for qG (0,00), 

EA«(i) - i'/(l+>/(A'log6)(c-l))-(l/log6)(A/p)(o''-l)^ 

and the Renyi function is given by 

T{q) =q(l + -^(e - 1)) - ^^-(6"^ - 1) - 1. 

Moreover, 

V.rAit)> r r (eV^°-"-' ^ 1) d.d. > 2t^ ' "'^^^ + ' 7 ^'^^^^^ + . 



Jo M 



1 /A' 



Proof. We have (5 = 1, and 

TTfc < const X 

fc! \/i 

Thus, 



5g :=^7rfc(c«'^'(Afc+Mfc) + fc) 

k=0 

V tl^'^iVM y ^(fc-l)!VM ; ^^(fc-l)!V/i 



< CX) 



for all positive A,/i,g. The condition for the inequalities (2.7) of Theorem 1 to hold then 
follows from Proposition 1. Now we consider the correlation decay and show that the 
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condition (2.5) of Theorem 1 holds. To see this, we consider the mother process 

A(t) = G(X(t)), G(u)=expju--(e-l) 

as a nonhnear transformation of the Markov process X{t),tE [0,1] which has marginal 
distribution 

p, = P{Xit) = j) = '"'^"j,^^^^' , J = 0, 1, 2, ... . 
Let C„(j; ^), 71 = 0, 1, 2, . . . , be Charlier polynomials, defined via the generating function 



(see, e.g., Chihara [7] and Schoutcns [29]). These polynomials form a complete system 
of orthogonal polynomials in the Hilbert space ^2(N, Poi(-)). In this case, the recursive 



for the Charlier polynomials: 



equation (3.4) for birth-death polynomials Qn{x) =C'n(^,^) is the recursive equation 



= nC„_i [ X,-] + ( X 71 ) C„ ( X, — ] + —Cn+l [ X, — ], 71 > 0, 



where 



^J.J \ Ai J \ fij fi 



Co{x,-]^l, C-i(x,- 



Then, the following expansion of the bivariatc distribution holds: 

P{X{t) = j,X{s) = fc) = p,p, e-'^-l*-!^ (j; ^) C„ (k; ^) 



n=0 

Note that 



G(7/)e/2( N,Poi( - 



since 



y g2(,-(A/M)(c-i))g-A/M (VmP ^ exp|--(c -l A< oo. 



We also note that G{u) has the expansion 
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au^Y^Ck (n; -] c-a/m (W:,„-(a/;.)(c-i) 1 



n=0 



Taking into account the orthogonality 



— di ^ k\ 



where dj. is the Kronecker delta function, we obtain 



Thus, 



and 



i?A(i, s) = Cov(A(t), A(s)) = ^c-^>l*l 

fc=i 



e-N^<i?A(r)<e--M5:^. 
'^i fe>i '^fc 

This completes the proof. □ 



4.2. Log-Pascal scenario 

B". Consider a mother process of the form 

where A'(t), t > is a stationary birth-death process with marginal Pascal- type distri- 
bution Pas(/3, ^) : 

Pj=P(X(0) = (m-A).7)= (l-^y^(^y, j = 0,l,2,..., 0<A<Ai, /3 > 0, 
and rates 

A„ = (n + /3)A, /.t„ = /in. 
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In the above expression for pj, we have used the Pochhammer symbol '■= P ■ (/3 + 
+ l),j> l,(/3)o = l. We note that 

-fG)' -M.^ 

Therefore the sum Sq in Proposition 1 is finite when q < log( j). 
Note that 

EX{t) = {^l- \)(3- / f 1 - -Y VarX(t) = - \f P- /{l-- 
Rx{t) = VarX(t)rx(t), rxit) = c"^*, t > 0. 
Under condition B", wc obtain the following moment generating function: 

M(C) =EeMCXm . < C < - log^/(. - A). 

We can formulate the following: 
Theorem 3. Suppose that condition B" holds. Then, for any 

1 - X/fi 



b > e"^'^^' 
the stochastic processes 



l-(l-A/^)e2(M-A) 



An(t)= / \\^^'\sV)ds, te[0,l], 



converge in Li to the stochastic process t G [0, 1] as 77, ^ oo such that, if A{1) G Lq 

for q e (0, min{log - log - A)}), 

w/iere 



log 6 

and i/ie Renyi function is given by 
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Moreover, 

rt rt 



A(l-c-''*) + l-e-^*(Ai< + l) 



Proof. We have again 6=1. In view oi q < log(^) and 



A' 

Sa < const X ki—c'^ 

the condition for the inequahties (2.7) of Theorem 1 to hold then follows from Proposition 
1. We now show that the condition (2.5) of Theorem 1 holds. For this purpose, we consider 
the mother process 

K{t) = G{X{t)), G(u)=exp{u-cx} 

as a nonlinear transformation of the Markov process X{t),t e [0, 1], which has marginal 
distribution pj ~ Pas(/3, ^) and the following bivariate distribution expansion: 

P{X{t) = J, X{s) = k)= p,p, V ^— A/„ /J, _ M„ fc; (3, - 



where Mn{j; /3, ^), n = 0, 1, . . . are Meixner polynomials, defined via the generating func- 
tion 



f:(^)„M„(,;/3,^)J = (l-.)--^(l 



-p , 0<z<- 



(see, e.g., Chihara [7] and Schoutens [29]). The recursive equation for the birth-death 
polynomial (3.4) becomes in this case the recursive equation for the Meixner polynomials: 

-jM„ (^j; 13, = bnMn+i (^j; 13, + 7„A/„ (^j; (3, + c„M„_i (^j;f3, j-^ n > 0, 



where 



and 



M_i( j;/3,^^ =0, Mo(^];l3,- ) = 1, 



A, ^, 1 n + {X/ii){n + P) n 

'---,^^ + ^^1^,' 1 A/, ' ^" = T^- 
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We note that 

G{u)el2 (N,Pas(^f3,^ 

since 

j=o V " 

We also note that G{u) has the expansion 



oo. 



°° / A\ 1 



/c=0 



P 



Taking into account the orthogonahty 

3'- VAi 



we obtain 



n,TO>l, 



A(t)=^afcA4(X(t);/3,-)^, 

oo 2 

i?A(i, s) = Cov(A(t), A(,s)) = ^ Je-'=^l*l 

fc=i "'fe 



and 



2 

e-Hj<i?A(r)<c-M;^ 

"'i k>i ' 



2 



This completes the proof. □ 

Remark 4- For A > /i, the birth-death polynomials are again given in terms of the 
Meixner polynomials: 
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and the stationary Markov process has the stationary Pascal-type distribution of the 
form 

p, = p{x{o) = (A - +0))=(^i^f^y^-^(^!^y, J = 0, 1, 2, ... , 

so that the conditions of Propostion 1 are satisfied for < g < log -. Thus, a scenario of 
the log-Pascal type can be produced for the geometric process X{t),t>0 by using again 
the orthogonality of the Meixner polynomials. We omit the details. 

For A = /i, the birth-death polynomials are given in terms of Lagucrrc polynomials 
L(f-i)(f),n = 0,l,2,..., as 

These polynomials are orthogonal with respect to the gamma density 

/(.) = 3^.'^- .>o. 

The corresponding log-gamma scenario can be produced, in principle, similarly to that 
given in Anh, Leonenko and Shieh [[3], [4]], but a proper interpretation is lacking in this 
case. 



4.3. Log-binomial scenario 

B'". Consider a mother process of the form 

A{t) = cMX{t)-cx}, cx=N\ogip{e-l) + l), (4.4) 

where X{t), t > is a stationary birth-death process with finite state space 5* = 
{0, 1, . . . ,N}, marginal binomial distribution Jiin{N,p): 

p, = P{X{0)^j) = (^^^ {f)p>{l-pf-\ J = 0,1,2, . . .,7V, <p< 1, 

and rates 

X„^{N-n)p, /i„=n(l-p), 0<n<iV. 

Note that 



EX{t)^Np, VarX(t) = 7Vp(l-p), 
Rx(t) = Y&vX{t)rx{t), rxit) = e-\t>0. 
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Under condition B'", wc obtain the following moment generating function: 
M(C) =Eexp{CX{t)} = (pe« + l-p)^, (eR. 
We can formulate the following: 

Theorem 4. Suppose that condition B'" holds. Then, for any 

(p(e^-l) + l)^ 
(p(e-l) + l)2^' 

the stochastic processes 

An{t)^ / r[A(J')(s6^')ds, te[0,l] 

converge in Li to the stochastic process A(t), t G [0, 1] as n oo such that, if A{1) € Lq 
for qe (0,oo), 

w/iere 

^('^) = '^(^ + g6)-l^^°^-^^(^'-^) + ^)' 
and i/ie Renyi function is given by 

^^^^-i'^^)-^b'^^^^^^'-'^^'^-'- 

Moreover, 

Y&rA{t)> [ [ exp(VarX(i)e-^l"-''l - l)dudw 
Jq Jo 

> 2tNp{l~p){{l - e-*) + 1 - c-\t + 1)). 

Proof. The condition for the inequalities (2.7) of Theorem 1 to hold follows from Propo- 
sition 1 since this is a process with a finite number of states, and thus Sq is finite for all 
q. Regarding the condition (2.5) of Theorem 1, we consider the mother process 

A{t)=G{X{t)), G(ii)=exp{ii-cjf} 

as a nonlinear transformation of the Markov process X{t),t G [0, 1], which has marginal 
distribution pj ~ Bin(iV,p) and the following bivariate distribution expansion: 

P(X(t)=j,X(s) = fc)=p,pfc^c-"l*-^li^„(j;7V,p)i^„(fc;iV,p) — , 

n—0 
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,2 _ {-irn\i{i-p)/py 



di = 



{-N)n 



where K„{j; N.,p),n = 0,1., ... , are Krawtchouk polynomials (see, e.g., Chihara [7] and 
Schoutens [29]). The recursive equation for the birth-death polynomial (3.4) becomes in 
this case the recursive equation for the Krawtchouk polynomials: 

-jKnU; N,p) = 6„if„+i(j; N,p) + 7„if„(j; N,p) + c„if„_i(j; N,p), n > 0, 
X_i(j;iV,p) = 0, Xo(j;iV,p) = l, 

and 

bn^p{N-n), 7„ = -(p(7V - n) + n(l -p)), c„=n(l-p). 

Note that 

G(M)eZ2(N,Bin(iV,p)), 

since 

N / at\ 

^e2(J-^-M ^ )p^(l-p)^-^ <oo. 

3=0 ^ 

We also note that G{u) has the expansion 

oo ^ 

G(M) = ^afeiC„(j;7V,p) — , 

k=0 



ak^Y. ^«(^'; N,p)r\pi{l- p)^-^e 



n=0 
N 



k = « j=0 

Taking into account the orthogonality 

TV 



j=o 

we obtain 



J2KnU;N,p)K^{r,N,p) (^^^ p^{l-p)^-^ = Sl^d: 



n, m > 1, 



RA{t,s) = Cov{A{t),A{s))^Y." 



d^ 

k=l "-k 
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and 



<i?A(T)<C-|^l^ 



A;>1 



4' 



This completes the proof. 



□ 



4.4. Log-hypergeometric scenario 

This subsection considers a birth-death process with quadratic rates. The hypergeometric 
function is meant to have the power series 

MaAc,z) = j:^-^^ = l + -z+ ^^^^^^ - + ..., 

where z is a complex variable and a, 6, c are parameters that can take arbitrary real 
or complex values provided that c ^ 0, — 1, — 2, . . . (see Abramowitz and Stegun [1] for 
details). 

B"". Consider a mother process of the form 

A(t)=exp{X(t)-cx}, cx=log-j^2Fi{~N,-g-h-N-e), (4.5) 

where X{t)^ t > is a stationary birth-death process with finite state space S = 
{0, 1, . . . , iV}, marginal hypergeometric distribution Hyp(A'^, h) : 

p,=P(X(0=j)-^4^' J = 0'---'^' 9.h>N, 

and rates 

\n = {N -n){g-n), fi„ ^ 7i{h - {N - n)) , n = 0,l,...,N. 

Note that 

EXit) ^ Var A(t) = ^i^/ia + + m^J + h - N) ^ 

g + h a + b — 1 

Rx[t) = Var A(i)rx(0, rx[t) = e"*, t > 0. 
Under condition B"", we obtain the following moment generating function: 

M(C) =Eexp{CX{t)} = i|l-2Fi(-A^, -g; h-N; e^), C € R. 
We can formulate the following: 
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Theorem 5. Suppose that condition B"" holds. Then, for any 

2Fii-N,-g;h-N;e^) 
{%)/&")[2F^{-N,^g;h-N;o^W' 

the stochastic processes 

converge in L2 to the stochastic process A(t), t G [0, 1] as 00 such that, if A{1) e Lg 
for qe (0,00), 

where 

\ ^ 1 r TP r AT „. J, i\T.„Q\\ 1 (n) 



^ii) = 91 + 7^-^ log(2^i i-N, -g; h-N;c'^))- log 



\ogbj \ogb "^^^ ' " a^')^ 

anii i/ie Renyi function is given by 



Moreover, 

Yai A{t) > f [ exp(VarX(i)e-l"-''l - l)ditdu 
^0 "'0 

^ jV(g/g + h)(l- g/{g + /i))(.9 + - TV) 
g + h-l 

(1 _ c-((g+/'+i)/2)t ^ 1 _ c-((g+/»+i)/2)t(((g ^ ^ ^ i)/2)t + 1) 
((5 + /i+l)/2)2 ■ 

Proof. The condition for the inequaUties (2.7) of Theorem 1 to hold follows from Propo- 
sition 1 since this is a process with a finite number of states, and thus Sq is finite for all 
q. In order to show that the condition (2.5) of Theorem 1 holds, we consider the mother 
process 

K{t) = G{X{t)), G{u)=eyip{u-cx} 

as a nonlinear transformation of the Markov process X{t),t& [0, 1], which has marginal 
distribution pj ~ Hyp((7, /i, TV) and the following bivariate distribution expansion: 

00 

P{X{t) = 3,X{s) = fc) = p,pk e-^(")l*-^li?,(A(n); g, h, N, )R,{\{n);g, h, N)p^, 
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A(n) = 71(71 — g — h — 1), 

_ {%-'-'){-NUm{-gU2n - g - h - 1) 
^" (-l)"?i!(-/i)„(n-.g-/i-l)jv+i ' 

where Ri are the dual Hahn polynomials defined by 

i?,(A(n)) = R,{X{n)-g, h, N) = Q^{i, g, h, N) = Q„(i), 

where Qnih g, h, N),n = 0, 1, 2, ... , arc Hahn polynomials defined by the following recur- 
rent relation: 

-xQn{x) ^ hQn+l{x) - (bn + Cn)Q„{x) + C„(3„_ i (.t) , Q-l{x) = 0, Qq{x) 1, 

where 

_ {n + g + h+l){n + g + l){N-n) _ n{n + h){n + g + h + N + 1) 

"~ (2n + 5 + /i + 2)(2n + .g + /i + l) ' ~ {2n + g + h){2n + g + h + 1) 

(see, e.g., Chihara [7] and Schoutcns [29]). Note that 

G{u)el2{^,Ryp{g,h,N)). 

The rest of the proof is similar to the proofs of the previous theorems using the orthog- 
onality of Hahn polynomials and expansion of the function G{X{t)) =Q^{'t)-'^x y^to a 
scries of Hahn polynomials with coefficients a„, n = 0, 1, 2, . . . , but in this case 

max A(?i) = max{7i(?i — g — h — 1)} 

n n 

is achieved for 

_ f g + h + 1 

and 

.2 °o „2 



'^1 k=i 
for dk ~ \/ JJfk- This completes the proof. □ 
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